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Abstract 

The goal of this paper is to sketch a broader outline of the mathematical struc- 
tures present in the Nonlinear Maxwell Theory in continuation of work previously 



presented in [11], \12\ and [13|. In particular, I display new types of both dynamic 
and static solutions of the Nonlinear Maxwell Equations (NM). I point out how 
the resulting theory ties to the Quantum Mechanics of Correlated Electrons inas- 
much as it provides a mesoscopic description of phenomena like nonresistive charge 
transport, static magnetic flux tubes, and charge stripes in a way consistent with 
both the phenomenology and the microscopic principles. In addition, I point at a 
bunch of geometric structures intrinsic for the theory. On one hand, the presence 
of these structures indicates that the equations at hand can be used as 'probing 
tools' for purely geometric exploration of low-dimensional manifolds. On the other 
hand, global aspects of these structures are in my view prerequisite to incorporat- 
ing (quantum) informational features of Correlated Electron Systems within the 
framework of the Nonlinear Maxwell Theory. 



1 Introduction 

The general goal of this paper is to examine broader ramifications of the Nonlinear 



Maxwell Equations (NM) as introduced by me in 1992/93 and further developed in JTT 



T2 1 , | I3|| . To this end, I first point out that the theory is considerably richer than that 



of the classical linear Electromagnetism. In particular, I describe here several distinct 
types of both static and dynamic solutions on a spacetime of the form M 3 x R. On 
the technical side, I have essentially avoided heavier analysis as the solutions are either 
obtained by means of elementary calculation, or are otherwise based on deeper analytic 
work described in [|13j. One should be aware that the possibilities opening in consequence 
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of the introduction of these new structures have not been fully exploited in this paper, 
thus postponing many potential developments into the future. 

More precisely, in the 'dynamic' part of the paper I display a solution in the form of 
a charge-carrying electromagnetic wave. It is a soliton type wave that transports charge 
with constant speed and without resistance. In addition, one notes existence of a specific 
to dimension four nonlinear Fourier type transform — an interesting structure whose role 
within the theory is twofold. On one hand, it can be used to find and analyze new solutions 
of the Nonlinear Maxwell Equations. On the other hand, the transform defines an exotic 
duality — a (quadratic) generalization of the (linear) Hodge duality. Consequences of this 
new duality for the four-geometry will be exploited in the future. 

The second set of results in this paper is focused around the question of existence and 
properties of static solutions. To this end, I first examine the situation on the Euclidean 
three-space. In particular, one takes note of the occurrence of global structures in the 
form of magnetic flux tubes as well as the so-called charge stripes. It is interesting from 
the point of view of geometry that these objects exist in general only on three-manifolds 
whose fundamental group is not finite. This is tied to the geometric fact that the nonlinear 
gauge theory at hand induces an additional structure on M 3 — namely a taut codimension 
one foliation. These global aspects of static solutions prompt an assumption of topological 
point of view. Accordingly, I sketch the possibility of constructing 'nonlinear cohomology' 
that would account for a sort of 'flux tube' invariant of a three manifold. The discussion 
here is based on two particular examples that I feel provide an optimal illustration of the 
underlying concept. 

The Nonlinear Maxwell Equations, cf. (0-0) below, involve a vector potential that 
encodes the electric and magnetic fields in the usual way as well as an additional scalar 
/. The function / contains information, extractable in a certain simple canonical way, 
about the local value of the filling factor (also known as the filling fraction). (The filling 
factor is defined as the number of quanta of the magnetic field per electron charge in the 
first Landau level. It is then natural and effective to think of the electrons as forming 
in conjunction with the corresponding magnetic flux quanta composite particles — either 
bosons or fermions, or Laughlin particles depending on the actual value of the filling 
factor.) It is thus postulated that the filling fraction — typically an input of a microscopic 
theory that is always assumed constant microlocally — is allowed to slowly vary in the 
coarser scale. In fact, it was shown in [O] that NM predict occurrence of phase changes 
that lead to formation of vortices in /, and a fortiori in the magnetic field. This picture 
conforms with the well known analogy between the Quantum Hall Effects and the High- 
T c Superconductivity. An inquisitive reader might now point at the following seeming 
conundrum. The physical interpretation of / as a filling factor requires the presence of 
two-dimensional geometric structures that endow us with a possibility of including the 
lowest Landau level in the basic dictionary. Thus, it may appear a priori puzzling, how we 
are going to retain this interpretation of / in three or four spatial dimensions? The answer 
is provided by the intrinsic structure of the NM themselves. On one hand, it is shown 
below that the filling factor variable may be completely factored out of the equations when 
viewed in the complete four dimensions of spacetime. Needless to say, if one attempted 
to analyze such /-free form in two dimensions the / variable would reemerge without 
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change as it is there encoded in the magnetic field B — bj f, b — const. On the other 
hand, one notes that a remnant or a generalization of the filling factor interpretation 
carries over to three dimensions. Namely, the NM in three dimensions imply existence 
of a codimension one foliation of the three-space associated with the static solutions. 
Moreover, one notes here that the NM do not a priori introduce any restrictions as to 
the type of the resulting foliation — in fact any regular foliation and even foliations with 
singularities introduced by degenerating leaves are admitted by the equations. However, 
as already mentioned above the existence of solutions of a special type, namely the flux- 
tube type, implies geometrical restrictions on the foliation and topological restrictions on 
the three-manifold. Indeed, in this case foliation must be taut. It also seems reasonable 
to expect that the composite particle interpretation remains valid in this setting and the 
number of participating electrons in each leaf is again determined by /, virtually leading 
to the notion of an effective Landau level. 

In the last words of this section I would like to admit that, the subject matter at hand 
being both new and inherently interdisciplinary as well as by way of my own background 
and limitations, it is not always easy to pick the optimal terminology. Realizing I will 
unavoidably fail to satisfy in this respect one group of readers or another, I can only ask 
the readership to be as tolerant as they can afford and hope that in the end substance 
will triumph over form. 



2 Nonlinear Maxwell Equations in Spacetime 

In what follows, in order to get around rather tedious algebra while not compromising 
our understanding of what is essentially involved, I present a shortcut style exposition of 
the necessary calculations. I believe that readers who are well familiar with differential 
geometry will find it easy to reinterpret this calculation in its natural invariant setting, 
while those who are less familiar with the abstract setting may in fact appreciate its 
absence here. 

Consider the following system of equations — the Nonlinear Maxwell Equations (NM) 
in the form in which they have appeared in my previous papers. 

dF A = (1) 

5(fF A ) = (2) 

nf + a\F A \ 2 f = uf. (3) 

where / is a real valued function and A is the electromagnetic vector potential, so that 
the corresponding electromagnetic field is Fa = dA. Here, a > is a physical constant 
with unit Tesla _2 m~ 2 ; I will not discuss the precise physical interpretation of a 2 in this 
paper. Further, d is the exterior derivative and 5 = *d* its adjoint. Here, it is assumed 
that the Hodge star * and the D'Alembertian □ are induced by the Lorentzian metric 
tensor on a spacetime of 3 + 1 dimensions. Let me point out that assuming / = const 
and dragging it to zero one recovers the classical Maxwell equations. In this sense, all 
phenomena of the classical electromagnetism are included in the present model. 
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Now the goal is to better understand the essential ingredients of the NM in terms of the 
classical field variables. To this end, let us say the spacetime is in fact the flat Minkowski 
space (with the speed of light 1) so that in particular one can identify coefficients of the 
electric field E and the magnetic field B with the coefficients of the curvature tensor Fa 
by the formula 

F A = Bidy Adz + B 2 dz A dx + B 3 dx A dy + E\dx A dt + E 2 dy A dt + E 3 dz A dt. (4) 

For the sake of our discussion below it is good to keep in mind the well known fact that 
the components of Fa are not Lorentz invariant. This property leads one to the derivation 
of the Lorentz force, so that the latter one is logically independent of the particular form 
of a gauge theory formulated in terms of Fa- In other words, the Lorentz force remains 
unchanged and valid as one attempts to modify the field equations. With this understood, 
let us continue the discussion of equations 

It would be rather straightforward to rewrite equations ([lflf) in the anticipated Maxwellian 
form by following the usual procedures for translating (|2|) into the Ampere and Gauss laws 
after having replaced E by fE and B by fB. In fact, this would lead to an ad hoc inter- 
pretation of / as a material constant — a route taken in our older, one might say naive, 
However, this form of the system offers little insight as to the more essential 
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paper 

implications of the NM, and one needs to find a less obvious reformulation. 
One notes that equation (0) may be equivalently written in the form 

f5F A = F A (V 3+1 f,.) 



(5) 



where V3+1 stands for the gradient in spacetime. For a reason that will become clear 
later, one identifies Fa with a skew-symmetric matrix in a standard way 
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It is important to note that by the miraculous property of skew-symmetric matrices in 
four dimensions, det F = E ■ B and 
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(I emphasize that F is not the matrix corresponding to the Hodge-dual of Fa in the given 
metric with signature (+ + H — ).) On the other hand, representing both the 1-forms and 
vectors as columns so that in particular 



3+1 
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and 5 F A 



Ei jt + B 2z — B 3jV 
E 2 ,t + B 3jX — B 1>z 
E 3) t + Bx tV — B 2)X 
E% >x + E 2 , y + E 3>z 
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one checks directly that 

^(Vg+l/, ■) = FV 3+1 f. 
This enables us to rewrite equation @ in the form 

F- 1 5F A = V 3+ iln(/). (6) 

It is perhaps worthwhile to realize that in this context F is a fiberwise- linear mapping from 
the tangent bundle to the cotangent bundle. Here one assumes / > a.e. This conforms 
with the principle that one will be consistently looking for strong solutions so that in 
particular / may always be replaced with |/| in ([IHI). Next one recalls that on one hand 
the first part of the NM (|IJ) is identical with the analogous part of the classical Maxwell 
equations and it encodes the Faraday's law of magnetic induction and the fact that there 
are no spatially extensive magnetic charges. This gives the first four (scalar) equations 
below, namely (0) and @. On the other hand, by direct multiplication and regrouping in 
(^), one obtains four further scalar equations that happen to radically modify the Ampere 
law. Written in the familiar three-space vector notation, the NM assume the form 

f)B 

V ■ B = (8) 
(— - Vxfl) xE + (V-E)B = -(E-B)V\nf (9) 



dt 



;^-Vx5).B = -(i.B^ln/ (10) 
(^-A)f + a(\B\ 2 -\E\")f = uf, (11) 



provided E ■ B ^ 0. In fact, also the case when E ■ B — is worth our attention and 
will be discussed below. As much as one should avoid indulging in formal manipulations 
of formulas, here the advantage of having the equations rewritten in several equivalent 
forms is that they all lead to the discovery of new types of solutions, the existence of 
which would be otherwise obscured by notation. This will become more evident in the 
following sections. 

As already mentioned in the introduction, I postulate the following physical interpre- 
tation: / is the spatially varying filling factor — a notion central to the modern composite- 
particle theories. In fact, the canonical microscopic-theory interpretation of the filling 
factor is valid in two spatial dimensions only, in which case it signifies the ratio of the 
number of quanta of the ambient magnetic field to the number of electrons in the first 



Landau level, cf. 0, fT7|. Moreover, the microscopic theory offers no hints as to the 



existence and relevance of an analogous notion in the three-space. A description of the 
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interaction of the electromagnetic field with fermions in the first Landau level provided 
by the equations above is valid in the mesoscopic scale. Here, as one 'zooms out' from the 
microscopic scale, the filling factor is neither a rational number nor is it a constant any- 
more. In fact, as it has been communicated in previous papers the spatially varying filling 



factor may assume the form of a vortex lattice, cf . |T3[ . For the time being, this point of 
view is validated by the well known analogy between the Quantum Hall Effect and the 
High Temperature Superconductivity and it awaits experimental confirmation. Moreover, 
the NM extend the notion of the filling factor to three spatial dimensions. However, as 
we will see below, the presence of the filling factor introduces an especially interesting 
modification of the laws of Electromagnetism only if the three-space comes equipped with 
a codimension one foliation. This latter fact makes it possible to talk about Landau levels 
in a certain sense, anyhow. Finally, / can be completely eliminated from the NM in 3 + 1 
dimensions. (In general, this requires that the first cohomology group of the spacetime 
vanishes.) In that case the NM can be written in the /-free form 

dF A = 

rf(F" 1 5F 4 ) # = (12) 

5 (F- l 5F^f - jF- 1 ^ 2 + a\F A \ 2 -u = 0, (13) 

where # is the isomorphism of the tangent and the cotangent bundles given by the metric. 
Indeed, under the assumption of vanishing first de-Rham cohomology, equation (^|) is 
equivalent to its integrability condition fll2|). Moreover, since the last scalar equation of 
the system can be written in terms of d In / in the form 

Win/- |rfln/| 2 + a|F A | 2 -z/ = 0, 

equation (^) also implies (|T3"D. Computation of the symbol shows that the system obtained 
in this way is non-hyperbolic — in fact its degeneracy is of higher order. Thus, this form 
of the NM appears impractical for any mathematical work, and an introduction of the 
dimensionless scalar / is necessary also from the point of view of analysis. Nevertheless, 
as indicated in the Introduction and the discussion above, physical implications of the 
existence of an /-free form of the NM are important. 



3 Geometry Behind the Equations 

The geometrical arena of the Maxwell equations consists of a spacetime, say N, and a 
principal £/(l)-bundle, say P, stack up above N. In addition, it seems any description 
of the interaction of the electromagnetic field with fermions requires, at least within this 
framework, a principal connection, i.e. a smooth (at least a.e.) distribution of horizontal 
planes that is invariant with respect to the circle action. This distribution can be written 
as ker A = ker fA for / ^ 0. In addition, if U(l) is to remain the elemental symmetry 
group of Electromagnetism, then / must be constant along the fibers so that it effectively 
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descends to a function on N. In particular, within this dictionary one can construct a 
Kaluza-Klein metric on P, which is given by 



Li A {X, Y) = g{ir*X, tt*Y) + aA(X)A{Y), 



Tesla 2 meter 2 



where the unit of a > must be 
[meter] and the unit of = dA is to remain, say, 
Laplace-Beltrami operator on forms is then A^ A = 
condition 

A A (fA) = vjA 
is equivalent to the system of equations (|l]-2D, cf. HTT| . 



if the unit of length on P is to be 
Tesla]. Let us say the corresponding 
= A^. Calculation shows that the 



4 Exotic Duality 



For the sake of discussion in this section, consider the NM on either a Lorentzian or a 
Riemannian four-manifold as the metric signature plays a secondary role. In particular, 
it is preferable to replace the D-notation with the A-notation. Assume for the sake of 
simplicity that the second cohomology group of the manifold is trivial. Omitting the 
constant a, write the system one more time in the form 

S(fdA) = (14) 

- Af+\dA\ 2 f = vf. (15) 
Since (|i"4|) implies d(f * dA) = 0, one has / * dA = dA so that 

dA = ±l/f*dA (16) 

and the new form A satisfies a dual system of equations 

6(jdA) = (17) 

-Af+\dA\ 2 j = uf. (18) 

This is a functional transform reminiscent of the Fourier or the Backlund transforms, 
notwithstanding the fact that all transforms are somewhat reminiscent of one another. 
In particular, the resulting dualistic perspective has the expected property that trivial 
solutions of one of the systems lead to more complex solutions of the dual system. To 
illustrate the idea, let me now present a few examples of dual solutions on R 4 with either 
the Euclidean or the Minkowski metric as specified in the discussion. 

Example 1. Let the metric be Euclidean and take dA = Edz A dt, E = const, and 
/ = f(x,y). Equation fli~4]) is automatically satisfied and fllBD assumes the form — f xx — 
fyy = i v ~ E 2 )f so that / = cos (k\X + k 2 y + a) for k\ + k\ = v — E 2 solves the problem. 
Now dA = ±f(x, y)dx A dy and it staisfies equations ([T7| - |T~8D . 
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Example 2. Departing for a while from the assumption of vanishing second cohomology, 
let us reinterpret the previous example on a four-torus assuming periodicity of coordinates 
(x, y, z, t) with period 2tc. Note that the first bundle is necessarily nontrivial as the 
cohomology class [dA] ^ 0. Let us allow the function / drop its dependence on y so that, 
say, / = cosx, provided the 'right choice' of v has been made. Now, dA = d(sinxdy) is 
an exact form so the second bundle is topologically trivial. 

Example 3. Consider dA = Bdx A dy and / = f(z,t) so that (|14D is satisfied. Let us 
now look at the metric with signature (+ + H — ) so that flT5| ) means f u — f zz — iy ~ B 2 )f. 
The general solution of this equation is a standing wave with variable amplitude. This 
pattern is inherited by dA = f(z, t)dz A dt (up to the sign again) which satisfies ( 

Example 4- Let us for a change begin on the other side and take, say, dA = edzAdt and 
/ = f(x,y). Again, the first equation ( |14|) is automatically satisfied while (p~5f) becomes 
—fxx — fyy + e 2 // = v f . As explained in |13| (see also remarks at the end of section |] 
below) apart from the trivial constant solution, this problem also has a solution in the 
form of a vortex lattice. In the latter case dA = f(x,y)dx A dy satisfies (P^-|T5|) and 
represents static magnetic flux tubes. 

I emphasize that only the vector potential A and the filling fraction variable / that 
appear in the first set of equations have physical interpretation. Reassuringly, the presence 
of a nontrivial / in examples 1 and 3 did not contribute anything unexpectedly strange 
to the constant electric and magnetic fields in these examples, while it 'introduced' flux 
tubes in example 4- Although one could consider similar interpretation of the transformed 
vector potential A, just as one can for any [/(l)-connection, I feel this is uncalled for and 
would probably be unjustifiable at this point. Nevertheless, the existence of the transform 
is a remarkable fact whose possible applications to four-manifolds will be explored more 
thoroughly in the future. In a way, this new duality is a generalization of the regular 
Hodge-star duality that may be compared to the projective generalization of the Euclidean 
reflection. This analogy may be justified in the following way. Projective duality is induced 
by a fixed quadratic form. What is the NM analog of that object? Introduce notation 
(p — Inf. A direct calculation shows that ([L4|-[15]) may be written in the form of a system 
of quadratic equations 

5dA + ~k(d(p A *dA) = (19) 

- A v - \dip\ 2 + \dA\ 2 - v = 0. (20) 

This form of the equation has one other advantage. Suppose one has found a solution 
(A, (p) of (|19| - |20|) . One can now use gauge invariance of the equations in the following way. 
Let x be a solution of the equation 

bd\ = -5 A. 

The existence of a solution x follows from the Fredholm alternative when the metric is 
positive definite, and it amounts to solving a linear wave equation in a Lorentzian metric. 
One can now replace A with A + dx (and denote the resulting form by A again). In the 
new gauge 5 A = 0, so that A in fact satisfies 

A A + *{d<p A *dA) = 0. (21) 
8 
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The system that consists of fl2"T| ) and fl2T)| ) is either quasilinear elliptic or hyperbolic, 
depending on the metric. Solving the latter system may not be helpful at all in finding 
solutions of the original ( pj| - |20|) , since one cannot guarantee that a solution satisfies the 
Lorentz gauge condition 6 A = 0. However, solutions of (|19|-p0|) a fortiori satisfy (^1|) and 
(p0|) so that in particular they will obey all a priori estimates on the solutions of, say, 
quasilinear hyperbolic systems. In particular, this point of view may justify the claim 
that the phenomena described in this paper shed some light on the complex nature of 
quasilinear systems of PDE of certain types in general. 



5 Charge Transport and Charge Stripes 



I will now take full advantage of the (|7)-|TT|) form of the NM. In analogy to the electromag- 
netic wave in vacuum, that one recalls is counted among the solutions of this system, one 
wants to look for a solution with E ■ B = 0. In the end I will check that the new solution 
of ([7|-pJ]) in fact satisfies (JH-^D which is not a priori guarantied. Make an Ansatz 



d d 
B = B 1 — + B 2 — , 
ox oy 



\ ox oy 



(22) 



where e, Bi and B 2 are a priori functions of (x, y, z, t) that are smooth a.e. and neither 
one of them vanishes identically. As an immediate consequence, one obtains that (0) and 
(§) are equivalent to 

Bt, t = {eB 1 ) >z 



B 



2,t 



(eB 1 ) >x + (eB 2 ) ;y = 
Bi >x + B 2t y = 







which implies 
On the other hand 



e, x B 1 + e„B 2 = 0. 



|) and ([10]) are equivalent to 

{B 2>x - B l>y )eB x = (eB^Bi - {eB x ) iV B x 

{B 2 , x - B 1 , y )eB 2 = (eB 2 ), x B 2 - {eB x ) jV B 2 
(-(eB 2 ) it + B 2>z ) B x + ((eBj),* - B 1>z ) B 2 = 0. 
and (p0|) imply that e is in fact constant 

e = ±1. 

Using (|23|) and fl2H|) again, one obtains 

B x = Bi(x,y,t + ez), B 2 = B 2 (x,y,t + ez). 



Equations (p3] 



(23) 
(24) 
(25) 
(26) 



(27) 

(28) 

(29) 
(30) 

(31) 
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In particular, B and E are not compactly supported. At this point, the only condition 
left a priori unfulfilled is the vanishing divergence condition. Thus, all equations ([23|-|3~0l) 

above are satisfied iff there is a function ip = ip(x,y,t + ez) such that 

Bi = -ijjy{x,y,t + ez), B 2 = ip x (x,y,t + ez). (32) 

Defining the electric and magnetic fields by (^) with e = ±1 , so that in particular 
\E\ = \B\, and choosing / that satisfies the linear wave equation (11), one obtains a 
solution of (|7|-[11]) . 

However, physical solutions must in addition satisfy the a priori more restrictive sys- 
tem (|l]-^D. Consider Fa as given in (|j). Equation (P is satisfied automatically since it is 
equivalent to (0-§|). On the other hand, (0) becomes 

(fB 2 ) >x - (fB 1 ) >y = (33) 

(feB 2 ) ;t -(fB 2 ) >z = (34) 
-(feB 1 ) )t +(fB 1 ) tX = (35) 
- {fBi),y + (fB 2 ), x = (36) 



Now, (j34|) and (|35D imply via that 

/ = f(x,y,t + ez). 

In particular, f jtt — f,zz = 0. Thus, (|l]-§) has been reduced to the following system of two 
equations: 

- f,xx ~ f,yy = Vf (37) 

(.A'.,)., • (/'••.«)., (38) 
The first equation above admits three types of classical solutions. Namely, 

( A(t + ez) In (x 2 + y 2 ) v = 

f=) A(t + ez) cos (hx + k 2 y + a(t + ez)) v = k{ + k\ (39) 
[ A(t + ez) exp (kix + k 2 y) v = —k\ — k\. 

Observe that each solution effectively depends on one harmonic variable in the (x, y)- 
domain — either, u = kix + k 2 y or u = lnr 2 = In (x 2 + y 2 ). Thus, equation (|38|) is satisfied 
if 

i) u = C(t + ez)/f(u,t + ez), 

for an arbitrary function C of one variable. Therefore, in view of (|32|) one obtains three 
types of solutions (redefining C) 

( C(t + ez)/(r 2 \nr 2 )[-y,x] 
[B u B 2 \ = \ C(t + ez) sec (fax + k 2 y + a) \-k 2 , k x \ (40) 
[ C(t + ez) exp (— k\x — k 2 y)[— k 2 , k\] 

in correspondence with ([39]). Since one is looking for strong solutions, one has the freedom 
to cut off pieces of the classical solutions (by restricting the domain) and to put them 
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back together. In this way, one obtains solutions that are either continuous or have 
jump discontinuities but may be guarantied to remain bounded. Last but not least, 
it is physically correct to interpret the divergence of the electric field as charge p and 
— ■§lE + V x B as the electric current. One checks that for solutions as above the (x,y)- 
component of current vanishes while the ^-component j is equal to — ep. More precisely, 
one obtains that piecewise 

f AC(t + ez)/(r 2 \n 2 r 2 ) 
ep = — j = < vC{t + ez) sec [k\x + k 2 y + a) tan [k\X + k 2 y + a) (41) 
[ —vC{t + ez) exp {—k\x — k 2 y) 

in correspondence with fl39f ) and (0). In addition to the piecewise smooth distribution of 
charge, one should include charge concentrated on singular surfaces where the electric field 
has jump discontinuities as indicated by the distributional derivative V ■ E. Therefore, 
charge is transported along the z-axis with the speed e = ±1 and without resistance as 
the vector of current is perpendicular to the electric field. Charge is mostly concentrated 
along charge stripes where the electric and magnetic fields have singularities. The net 
current depends on the particular choice of a (strong) solution. Of course, the theory 
does not tell us how to solve the practical problem of electronics — namely, how to create 
conditions for a particular function C = C(t + ez) , constant v and a desired mosaic of 
singularities to actually occur in a physical system. 



6 Static Solutions and Magnetic Flux Tubes 

The classical Maxwell equations admit static solutions of two types only: the uniform field 
solutions, and the unit charge or monolpole-type solutions, as well as superpositions of 
these fundamental types of solutions. As we will see below, the nonlinear theory encom- 
passes a larger realm including the magnetic-flux-tube type and the charge-stripe type 
solutions. These additional configurations require nonlinearity and cannot be superposed, 
which gives them more rigidity. In the next section we will see what can be said about the 
variety of such solutions, while in this section I will only display a single example of this 
type. Apart from the applicable goal, the idea is to present an example that possesses all 
the essential features of the general class of solutions yet the required calculation is free 
of more subtle geometric technicalities. 

Time-independent solutions of the NM posses physical interpretation only if they 
satisfy the equations in the classical sense almost everywhere. Assuming that all fields 
are independent of time (|7[-|ll"|) takes on the form 



V x E = (42) 

V • B = (43) 
(V x B) x E+ (V -E)B = -(E- B)V\nf (44) 

{V x B)- B = (45) 



11 



-Af + a(\B\ 2 -\E\ 2 )f = uf, (46) 

under the assuption that E ■ B ^ a.e. Adopt an Ansatz that the integral surfaces of the 
planes perpendicular to the field B are flat, say, 

d 

One easily checks that equations ( ^3|) and ( f4"5|) are satisfied. Assume in addition that the 
electric field is potential, i.e. 

E = S/ip(x, y, z), where ip z ^ a.e. 

so that (|42| ) is satisfied. Remembering notation = In /, one calculates directly that 

-> -> d d d 

(V x B) x E = -ip z b x — - ^ z b y — + (ip x b x + ip y b y )—, 

while 

(V-E)B = A#^, 

and 

(E ■ B)V<p. 

Thus, equation (|44| ) is equivalent to the following system of three equations 

ip z (b(p x + b x ) = 

i/) z (bp y + b y ) = 
b Aip - ip x b x - ipyby - bip z ip z = 
and since ip z ^ one obtains from the first two equations 

ip(x, y, z) = ipx(x, y) + <p 2 (z) and b = (3 exp (-y?i), 

while the third equation assumes the form 

Aip + Vip ■ V<p = (47) 

At this point the NM have been reduced to the system of just two scalar equations (|4*B] ) 
and (|47j). Denote fi = exp(pi and / 2 = exp<^ 2 and assume in addition 

ip = ip(z) 

so that 

^i'(z) = eexp(-</? 2 ) = v- 

J2 

It now follows from ( fiq) and fl4"7|) that the triplet 

h{x,y)oz j2{z)az 
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and 

f(x,y,z)=f 1 (x,y)f 2 (z) (49) 

is a solution of the NM if only f\ and f 2 satisfy a decoupled system of semi-linear elliptic 
equations 

-tf(*)-7^r = (50) 

- A/i(z, y) + ^ = i^/i(ar, y). (51) 

At this point, I would like to emphasize one more time that in a field theory one 
looks for strong solutions, i.e. solutions that satisfy equations in the classical sense almost 
everywhere. Typically, such solutions are smooth except for singularities supported on 
a union of closed submanifolds. Furthermore, geometrically invariant derivatives of the 
resulting fields in the distributional sense signify charges. With this understood, let us 
briefly turn attention to equation One wants to avoid holding the reader hostage 

to the formal analysis of this elementary equation which might be somewhat distracting. 
Thus, I have chosen to briefly describe the solutions qualitatively leaving aside technical 
details that can be easily reconstructed aside by the reader. First, one notes that if 
u 2 > then a solution is concave, while for v 2 < it will be convex for large values where 
/f > —e 2 /v 2 . Assuming formally that f 2 is a function of f 2 (piecewise), one reduces fl5"Cf ) 
to the first order equation 



Thus, there are essentially two types of positive solutions, depending on the actual values 
of constants c,e,u 2 . The first type includes solutions that assume value at a certain 
point Zq and increase monotonously to infinity as z — > oo as well as the symmetric 
solutions defined between — oo and some point, say zq again, where they reach 0. These 
solutions require v 2 < and they asymptotically look like exp (±(— i> 2 ) l l 2 z) One can use 
both branches in order to put together a strong solution that forms a cusp or a jump 
discontinuity at zq. The second type consists of solutions that are concave, rise to the 
highest peak at f 2 = m, when c — v 2 j\ — £ 2 ln/f = 0, and fall off to on both sides in 
finite time while being differentiable in-between. Selecting the constants and combining 
both types of solutions piecewise segment-by-segment one obtains strong solutions f 2 that 
in turn provide electric fields according to formula fl4~8[) . 

Since, with the exception of the trivial constant solution, there are no global smooth 
solutions, one concludes that either E is constant or there exist charge stripes located at 
planes z = const where f 2 (z) has singularities. The distributional derivative is in each 
case equal to the Dirac measure concentrated at z — const as above and scaled by the size 
of the jump, and classical derivatives on both sides of the singularity. Even in absence 
of a jump, the charge will switch from negative to positive thus forming what can be 
amenably called a charge-stripe. An example of this is shown in Fig.l. 

It is much more difficult to figure out solutions of the second equation. I refer the 



reader to JT3[ for a more thorough analysis, while here I will just briefly summarize my 
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previous findings. Solutions of equation Q5ip correspond to critical points of the functional 



y\vA\ 2 +{3 2 fHh: 



If! 

which is neither bounded below nor above, so that one is looking at the problem of ex- 
istence of local extrema. The equation always admits a trivial constant solution. But, as 



it is shown in fL3|] , it also possesses nontrivial vortex lattice solutions. More precisely, if 
(3 is larger than a certain critical value then there is a nonconstant doubly periodic func- 
tion / which satisfies the finite difference version of Q5ip everywhere except at a periodic 
lattice of isolated points, one point per each cell. In this way, a lattice of flux tubes, 
cf. Fig. 2, emerges as a solution of the NM. For the time being, the proof of this fact 
relies on finite-dimensionality essentially, and does not admit a direct generalization to 
the continuous-domain case. However, physical parameters, like / f 2 and (3, are asymp- 
totically independent of the density of discretization. Thus, I conjecture existence of the 
continuous domain solutions that satisfy the equation a.e. in the classical sense and re- 
tain the particular vortex morphology. Presently, the essential obstacle to proving this 



conjecture is lack of a regularity theory for the discrete vortex solutions. The proof in [|L3] 
is carried out in the (discretized) torus setting. One believes that vortex type solutions 
exist on any closed (orientable) surface. 



7 Topological Quantum Numbers 

Every gauge theory comes equipped with an associated set of topological invariants— 
usually characteristic classes of the bundles used to introduce the gauge field. Articles J|] , 
||, teach us how such topological invariants may be manifested in an electronic system 
as observable quantum numbers. The Nonlinear Maxwell Theory is naturally equipped 
with two kinds of topological invariants. On one hand, one has the first Chern class of the 
original [/(l)-bundle. Additionally, we will see below that in the case of static solutions 
the NM give us an additional set of invariants defined directly by the foliated structure 
of the underlying three-manifold. (In the discussion below, I generally assume for the 
sake of simplicity that M is a closed orientable manifold unless stated otherwise.) In 
this section I will make an effort only to identify rather than exploit to the fullest the 
geometric and topological ramifications of this nonlinear theory of Electromagnetism. To 
gain some initial impetus, let us be guided by the following question 

What are the necessary and sufficient conditions on a Riemannian three- 
manifold M for the NM to admit a separation of variables of the type 
seen in the previous section, i.e. for the equation to decouple so that 
its solutions will generate magnetic- flux-tube type solutions on M? 

A question of this type is typical in algebraic topology where one is asking about global 
obstructions to the presence of certain algebraic factorization properties of analytic ob- 
jects, like linear differential equations as it is the case for, say, the de-Rham cohomology 
groups. In our case, the equations are nonlinear, but the principle remains the same. The 
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importance of these questions for practical issues of Electromagnetism is twofold. First, 
one wants to know how big is the set of possible configurations — especially in the absence 
of the superposition principle. Secondly, I believe the topological invariants displayed 
below are directly on target in an effort to explain and describe the nature of certain rigid 
structures, like the Quantum Hall Effects, that physically occur in electronic systems. 

First, it needs to be emphasized that the static field equations I want to consider, 
i.e. the equations that descend from the four- dimensional spacetime via time-freezing 
coefficients, are distinct from the equations ([jJIJ) considered directly on a three manifold. 
Secondly, the equations fl42| - |4"6"D are only valid on a Euclidean space. The geometry behind 
these equations is easier to identify when they are rewritten in an invariant form that can 
be considered on any three-manifold in a coordinate independent setting. 

Fix a Riemannian metric on M with scalar product < ., . > extended to include 
measuring differential forms. Denote by B and E the forms dual to the magnetic and 
electric field vectors; recall notation (p = In / and put a = 1. The static NM assume the 
form 

dE = (52) 

SB = (53) 

★ (*dB A E) + (5E)B = - < E,B > dtp (54) 

dB A B = (55) 

A^+|c^| 2 + |£| 2 -|B| 2 + z/ = 0. (56) 

Equation (|55|) is the familiar Frobenious condition on integrability of the distribution of 
planes given by ker B. One always assumes B is nonsingular a.e. so that the distribution 
is a priori also defined a.e. For convenience, it is assumed throughout this section that 
the foliation determined by ker B is smooth. (It is quite clear that for the flux-tube type 
solutions the distribution extends through the singular points and is defined everywhere. 
At this stage, however, it is hard to make a formal argument to this effect, hence the a 
priori assumption.) The condition of smoothness implies that the three-manifold M must 
have vanishing Euler characteristic. In particular, singular foliations, some of which may 
be associated with other types of solutions of the NM, are excluded from the discussion 
below. 

It follows that there is a 1-form a, known as the Godbillon-Vey form, such that 

dB = a A B. 

This form is not defined uniquely. However, as is well known, d(a A da) = and the 
Godbillon-Vey (GV) cohomology class 

[a A da} H 3( M ) 

is uniquely defined. On a three manifold this class can be evaluated by integration result- 
ing in a Godbillon-Vey number 

Q — a A da. 
Jm 
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This invariant poses many interesting questions that have not been fully resolved by 
geometers yet. Below, I will justify two observations. First, the condition of existence 
of the magnetic flux-tube solutions imposes both local and global restrictions on the 
foliation. Second, magnetic flux-tube solutions exist in topologically nontrivial situations 
with nonzero GV-number Q. This is formally summarized in the two propositions that 
follow. They are far from the most general statements that can be anticipated in this 
direction, but are also nontrivial enough to suggest a conjecture regarding quantization 
of the GV-number that I will formulate following Proposition [I]. 

Consider a priori a foliation given by kerS locally. First, one introduces a local 
coordinate patch (x, y, z) such that the foliation is given by the (x, y)-planes and \dz\ = 1. 
In particular 

B = (3(x, y, z)dz. 

Let 7 = g(x, y, z)dx A dy denote the volume element on a leaf. One has *B = f3(x, y, 2)7. 
Equation (|53|) becomes d(j3(x,y, z)g(x,y, z)dx A dy) = 0. Thus, there is a function x — 
x{x, y) such that 

g{x,y,z) 

A calculation analogous to that in the previous section shows that the whole system 
(pM)B[) is reduced to 

7 =~<Px, In -7 r =-<p y (57) 

V 9{x 1 y,z)J x \ g{x,y,z)J y 

5E+ < E,dip>=0 (58) 



2 



Aip+ \dip\ 2 + \E\ 2 -[ , +v = 0. (59) 

\9{x,y,z)J 

Observe that in order to obtain factorization 

<p(x,y,z) = <fi(x,y) + ip 2 (z) (60) 
it is necessary and sufficient that 

9 = 9(x,y), (61) 

i.e. a priori dependence of g on z is dropped. If that holds, the equations fl58f ) and fl59| ) can 
be decoupled with an additional Ansatz E = e(z)dz. One also has that %jg = 6exp (— tpi) 
for a constant b and 

A x , y y?i + \difi\ 2 - 6 2 exp (-2</?i) + v = 0. (62) 

Conversely, if (32) and (|60D hold, then so must d6l|) and the mean curvature h of a leaf 
vanishes. Indeed, by definition 

h — S [t^tB] = — -k d(g(x, y)dx A dy) = 0. 



LB 



This implies 
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Proposition 1 For the existence of flux-tube type solutions — in the sense of existence of 
factorization and decoupling of equation (ffi&J — it is necessary that the foliation given 
by ker B be taut, i.e. the mean curvature of leaves must vanish. In particular it \(M) must 
be infinite. 

Proof. The first part has been shown above. The second part follows from a result of D. 



Sullivan that he deduced from the result of Novikov on the existence of a closed leaf 



that is a torus (cf. ||, and |TJ| for additional general material and references) . □ 



In particular, there are no flux-tube type solutions of the NM that would conform 
with the Reeb foliation |J. This is a practical issue since the Reeb foliation exists on 
a solid torus, so that in principle it might be observed experimentally which would be 
inconsistent with the theory at hand. This fact is also interesting for another reason. 



Namely, according to the celebrated theorem by Thurston in |L5[] each real number may 
be realized as the Godbillon-Vey number for a certain codimension one foliation on the 
three-sphere S* 3 . The known proof of this result uses the Reeb foliation in an essential 
way. I do not know if this fact is canonical, i.e. if the presence of the Reeb foliation 
is necessary for the result to hold, but if it turns out to be so then excluding the Reeb 
foliation from the game should result in a reduction of the range of the G-V number, 
possibly to a discrete subset of the real line. In such a case, the resulting set of the 
G-V numbers accompanying flux-tube type solutions of the NM would also be discrete. 
This is consistent with my expectation that these invariants must be related with (both 
the integer and the fractional) Quantum Hall Effects. Future research should bring a 
resolution of this problem. 



Another observation is that the factorization given by (60) and fl62|) does exist in 
topologically nontrivial situations. More precisely, I want to consider solutions of the 
NM on PSL(2, R) and its compact factors. These three-manifolds are equipped with 
interesting codimension one foliations known as the Roussarie foliation [|10|| . Let the Lie 
algebra SL(2, R) be given by 

[X,Y + ] = Y + 1 [X,Y~]=-Y~, [Y + ,Y~]=2X. 

Pick a metric on PSL(2, R) in which the corresponding left-invariant vector fields X, Y + , 
and Y~ are orthonormal and let fi, v + , and u~ be the corresponding dual 1-forms. One 
checks directly that 

dv~ = n A v~ . 

so that the distribution ker v~ is integrable and fi is the GV-form of the resulting fo- 
liation. In particular, one can introduce local coordinates (x, y, z) such that d x = X, 
d y = Y + , d z = Y~ . This foliation descends to compact factors of PSL(2,R) that can 
each be identified with T\M g — the total space of the unit tangent bundle of the hyper- 
bolic Riemann surface of genus g that depends on our choice of the co-compact subgroup 
acting on PSL(2,R) by isometries. Moreover, the GV-integrand /i A dfi is proportional 
to the natural volume form on the three-manifold. As a result of this, the corresponding 
GV-numbers 

Q=[ nAd/i = -2Vol(M q ) 

JTlMg 
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assume values in a discrete set. I want to look for solutions of the NM that satisfy the 
Ansatz 

B = (3is-. (63) 

In particular, (the Frobenious) equation (|55D is satisfied automatically. Moreover, since 
•kdB = (Y~j3)fi A u + A u~, equation (|53| ) implies 

(3 = (3{x,y). (64) 

As before, one checks that Q5"3] ) implies ( plf ) as well as X(p = —X In (3 and = 
— y + ln/5. In consequence, one again has fl6"0|) and assuming = e(z)^~ as before one 
obtains flBlq) . In consequence, the following holds true. 

Proposition 2 The Roussarie folitions on PSL(2,R) and its compact factors satisfy the 
factorization condition for the existence of magnetic flux-tube type solutions in the sense 
that the tangent distribution can be expressed as ker B a.e. and one can reduce the NM 
to the form (fg^-[g^j. 

In a similar way one can obtain factorization ( |5"D| ) and ( |B"2D for other foliations, like the 
natural foliation on say S 2 x S 1 . 



8 More on the Physical Framework of the NM 

It is natural to ask if the NM descend from a Lagrangian functional depending on the 
two variables A and /, say <&(A,f), via the Euler-Lagrange calculus of variations. The 
answer is negative as one can easily see considering that in general a gradient must pass 
the second derivative test: 

5 2 S 2 
$ = $ 

SASf 5f8A 

—a condition that cannot be satisfied by the expressions in (|l]-|D viewed as the gradient, 
say °f an unknown functional $. This suggests that the NM may constitute 

just a part of a broader theory that would encompass additional physical fields. In other 
words, the equations ([i]-||) would have to be coupled to some other equations via additional 
fields. In addition, such coupling would have to induce only a very small perturbation of 
the present picture that one believes is essentially accurate. Such possibilities may become 
more accessible in the future. Among other, perhaps related goals is that of deriving the 
NM equation directly from the microscopic principles. 

The well-known analogy between the Quantum Hall Effects and High-T c Supercon- 
ductivity suggests that there should exist vortex lattices involving the so-called filling 
factor (microlocally a constant scalar) that plays a major role in the description of Com- 
posite Particles. The NM describe exactly this type of a vortex-lattice. Simulation and 
theory show that this system conforms with the experimentally observed physical facts. 
It stretches the domain of applicability of the Maxwell theory to encompass phenom- 
ena such as the Magnetic Oscillations, Magnetic Vortices, Charge Stripes that occur in 
low-temperature electronic systems exposed to high magnetic fields. 
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There are other systems of PDE that admit vortex-lattice solutions and are conceptu- 
ally connected with Electromagnetism, like the well known Ginzburg-Landau equations 
valid within the framework of low T c type-II superconductivity, or the Chern-Simons ex- 
tension of these equations which, some researchers have suggested, may be more relevant 



to the Fractional Quantum Hall Effect and/or High-T c Superconductivity, cf. [IS]. The 
free variables of these equations are the so-called order parameter (a section of a com- 
plex line-bundle) and a [/(l)-principal connection, both of them containing topological 
information. In the case of NM, all the topological information is contained in one of the 
variables, i.e. the principal connection, while the other is a scalar function. An additional 
advantage of the NM is in that it remains meaningful in three-plus-one dimensions just 
as well as in the two-dimensional setting. I would also like to mention that recently other 
researchers have introduced nonlinear Maxwell equations of another type in the context 
of the Quantum Hall Effects, cf. ||. The NM theory presented in this and the preceding 
articles of mine is of a different nature. Finally, although this is far from my areas of 
expertise and the remark should be received as completely ad hoc, I would also like to 
mention that yet another context in which foliations come in touch with the Quantum 
Hall Effect is that of noncommutative geometry, cf . M . 

Let me conclude with a question that may suggest yet another point of view. Namely, 
is there a coalescence between the nonlinear PDEs (in the form of the NM) and the 
(Quantum) Information Theory? As it was pointed out, construction of error correcting 
codes may unavoidably require manipulating quantum information at the topological 
level. Anyhow, this is how I have understood the essential thought in M. Adopting this 
paradigm would strongly suggest that the effective language of quantum computation 
should be costructed at many levels, including that of the mesoscopic field theory in 
parallel with the language derrived from the basic principles as it is done now. Future 
research will likely better clarify these issues. 
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Fig.l An example of a strong solution of Q50|). / = is a positive function, the 

electric field is given by formula (0). The resulting charge distribution is obtained by 
evaluating V • E. (In general, V • E is understood in the distributional sense). Charge 
is concentrated along certain plains z = const. This is the basic appearance of charge 
stripes — intertwining concentrations of positive and negative charges. (One should com- 
pare this static picture with the description of moving charge stripes in section |^.) 



Strong solution f=f(z) 




The corresponding electric field along the z-axis 




The resulting charge distribution along the z-axis 
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Fig. 2 The luminance graph of / = fi{x,y) that solves ([H]). The corresponding magnetic 
field on the right is obtained via ([48]). 
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